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1 Introduction 



The measurements of the deep-inelastic scattering (DIS) structure function (SF) F2 in 
HERA H, Q have permitted the access to a very interesting kinematical range for testing 
the theoretical ideas on the behavior of quarks and gluons carrying a very low fraction 
of momentum of the proton, the so-called small x region. In this limit one expects 
that nonperturbative effects may give essential contributions. However, the resonable 
agreement between HERA data and the NLO approximation of perturbative QCD has 
been observed for > IGeV^ (see the recent review in [3) and, thus, perturbative QCD 
could describe the evolution of structure functions up to very low values, traditionally 
explained by soft processes. It is of fundamental importance to find out the kinematical 
region where the well-established perturbative QCD formalism can be safely applied at 
small X. 

The standard program to study the small x behavior of quarks and gluons is carried 
out by comparison of data with the numerical solution of the Dokshitzer-Gribov-Lipatov- 
Altarelli-Parisi (DGLAP) equations |§] by fitting the parameters of the x profile 
of partons at some initial Qq and the QCD energy scale A p[-|Tl|. However, if one is 
interested in analyzing exclusively the small x region, there is the alternative of doing 
a simpler analysis by using some of the existing analytical solutions of DGLAP in the 
small X limit This was done so in Ref. [l^, |l3l where it was pointed out that the 



HERA small x data can be interpreted in terms of the so called doubled asymptotic scaling 
phenomenon related to the asymptotic behaviour of the DGLAP evolution discovered in 
1^, many years ago. 

On the other hand, various groups have been able to fit the available data (essentially 
at large Q^) using a hard input at small x: A > 0. In some sense, it is not very 

surprising, because the modern HERA data cannot distinguish yet between the behavior 
of a steep input parton distribution and the quite steep dynamical evolution from a soft 
initial condition. Moreover, for the full set of anomalous dimensions (AD) obtained at 
X — ^ in Ref. based on BFKL, the results weakly depend on the form of the 

initial condition (see [^), preserving the hard ones and changing the soft ones. In the 
case of working with fixed order AD, the initial conditions are important when the data 
are considered in a wide range of and it is necessary to choose the form of the PD 
asymptotics at some Qq. In this work we use a soft initial condition in agreement with 
the experimental situation: low-Q^ data P^l E31 are well described at < OA GeV^ by 



Regge theory with Pomeron intercept ep = X + 1 = 1.08, closed to the standard {ep = 1) 
one. Moreover, HERA data with > 1 GeV^ are in good agreement with GRV 
predictions |lT| which support our aim to develop an analytical form for the parton 
densities at small x because, at least conceptually, it is very closed to the GRV approach. 

Thus, the purpose of this article is to obtain the small x asymptotic form of parton 
distributions (PD) in the framework of the DGLAP equation starting at some Qq with 
the fiat function: 

UQl) = {a = q,g), (1) 

where fa are the parton distributions multiplied by x and Aa are unknown parameters 

^ At small X there is another approach based on the Balitsky-Fadin-Kuraev-Lipatov (BFKL) equation 
||7|, whose application is out of the scope of this work. 
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that have to be determined from data. Through this work at small x we neglect the 
non- singlet quark component. 

The article is organized as follows. Sections 2 and 3 contain LO analyses: in Sect. 2 
we consider the simple but important case without quarks, while in Sect. 3 the quarks 
distributions are taken into account. The NLO analysis, which is the main result of this 
article, is performed in Sect. 4. In Sect. 5 we present the fits to experimental data and 
some discussions of the obtained results. In the Appendix we illustrate the method of 
replacing at small x the convolution of two functions by a simple product. The method is 
used in the present work for the correct incorporation of the non-singular part of parton 
distributions to our formulae. 



2 Leading order without quarks 

First of all, we consider the leading order (LO) approximation without quarks as a 
pedagogical example of the more cumbersome calculations below. This case is at the 
same time very simple and very closed to the real situation, because gluons give the basic 
contribution at small x. 

At the momentum space, the solution of the DGLAP equation in this case has the 
form 



Mg{n,Q') = M,(n,Qg)e-'^..W^ (2) 
where Mg{n, Q^) are the moments of the gluon distribution. 

The terms lfg{n) and /3o are respectively the LO coefficients of the gluon-gluon AD and 
the QCD /3-function. Through this work we use the short notation a(Q^) = a;s((5^)/(47r). 
At LO, s can be written in terms of the QCD scale A as: 



For any perturbatively calculable variable K{n), it is very convenient to separate the 
singular part when n ^ 1 (denoted by "i^") and the regular part (marked as "i^'"). Then, 
Eq. (0) can be represented by the form 



Mg{n,Q'^) = M3(n,Q2)e-49^Lo/(n-i)e-'i9.W«Lo^ (4) 

with ^gg = -8Ca and Ca = N for SU{N) group. 

Because we are only interested in the small x behaviour and the initial conditions are 
given by the soft (x-independent) functions in Eq. (|l]), we use permanently the variable 
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z = x/xo with values < 2; < 1 with some arbitrary Xq < 1 0, i.e. 

M,(n,g2) = f'dzz^-^Uz,Q^) 
Jo 

Finally, if one takes the soft boundary conditions given by Eq. (0), the coefficient in 
Eq. (0) becomes 

Main,Ql) = (5) 
n — I 

2.1 Classical double- logarithmic case 

As a first step, we consider the classical double-logarithmic case which corresponds to 
use only the AD singular part, i.e. dgg{n) = in Eq. (3). 

Then, expanding the second exponential in the r.h.s. of Eq. (3) 

k 

^""'^ ^ - ^'^^k\ (n- 1)^+1 



and using the Mellin transformation for (ln{l/ 



k 



dzz''-^{ln{\lz)f 



(n - 1)^+1 

we immediately obtain the well known double-logarithmic behavior |T9| : 

ff{z,Q^) = AgJ2--i-dggSLo) {Hl/z)f = AgI,{aLo). (6) 
fc=o y'^-i 



where Io{<^lo) is the modified Bessel function with argument^ cr^o = "^ydggSLo^n^z). 
2.2 The more general case 



For a regular kernel K(z) (see and Appedix), having Mellin transform 



K(n) = dzz^'-^kiz) 
Jo 



and the PD fa{z) in the form dln{l/ z)) we have the following equation 



K{z)®Uz) = K{l)Uz)+OUj^^^ 



(7) 



^ The correct incorporation of the regular part of the parton distributions can be done only for small 
values of z (see Eq. (^) . Thus, xq should be restricted to the range 0.1 < < 1 in order to keep the 
correctness of our formulae at x < 10'^. 

^Hereafter we follow the popular Ball-Forte variables ctlo and pLO- Below, they are generalized to be 
used beyond the LO approximation (cr and p). 
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The Eq. has been obtained for the nonsingular (at n — 1) functions K{n) and 
Regge-hke PD many years ago by the Madrid group ||25|. It has been expanded for 



arbitrary kernel K{z) in p4[ and it was aheady used in pGf and to extract the gluon 
distribution and the longitudinal structure function from F2 and dF2/ dln{Q'^) data. 



For arbitrary K{z) the ability of the method |24| has been checked numerically in Refs. 
|2§ and using MRS sets of PD. 

With the help of Eq. one can find the general solution for the LO gluon density 
without the influence of quarks 



fg{z,Q'^) = AgIo{aLo)e 



-dgg{l)SLO 



0{p 



LO) 



where 



Plo 



dggSLO _ O'LO 

\ ln{z) ~ 2ln{l/z) 



, 7^(1) = 22 + ^/ and dgg{l] 



1 + 



3/3o 



with / as the number of active quarks. 



3 Leading order (complete) 

At the momentum space, the solution of the DGLAP equation at LO has the form 



M+{n,Q^) + M-{n,Q^) and 



■d±s/{n~l) -d±{n)s 



(9) 



where0 



d±{n) 



2/3o ' 



1 + 



Adqg{n)dgq{n) 



{dgg{n) + dqq{n)) ± {dgg{n) - dgq{n 

Vi I dM-dgg{n) \ ^ 
^99^^') 2 I ^ d±{n) - d^{n) )' '''^ ^ d±{n) - d^{n) 



\ {dgg{n) - dqq{n)) 

dab{n) 



a ^ 6110) 



As the singular (when n ^ 1) part of the + component of the anomalous dimension is 
d+ = dgg = —4:Ca/Po while the — component does not exist {d- = 0), we consider below 
both cases separately. 



^We use a non-standard definition (see psf) of the projectors s^i^in), which is very convenient beyond 
LO (see Eq. (|2^)). The connection with the more usual definition a, a and e in ref. |2^, |3^ is given by: 
£<7<,(") = ain), e- (n) = a{n) and (n) = e(n) 
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3.1 The "+" component 

The analysis of the "+" component is practically identical to the case studied in section 
2. The only difference lies in the appearance of new terms B'^fj{n). If they are expanded 
in the vicinity of n = 1 in the form ^^^(n) = + (n — l)^^^, then for the terms 
multiplying Mb{n, Q^), we have the same results as in previous section: 

where the symbol — > denotes the inverse Mellin transformation. The values of a and p 
coincide with those defined in the previous section because (i+ = dgg. 

The terms that come with the additional factor (^—1) in front, lead to the following 
results 



k 

(rj - _A^^-d+SLo/{n-l) _~+ a 1 {-d+SLo) 



i.e. the additional factor (^ — 1) in momentum space leads to replacing the Bessel function 
Io{o-Lo) by pLoh{(TLo) in 2-space. 

Thus, we obtain that the term £^b(n)Mft(n, Q^) leads to the following contribution in 
z space: 

elM<yLo)+et,PLoh{(rLo)]A,e-~^+^^'^'^o + o{plo) (H) 



oo 



1 1 



Because the Bessel function Iv{ct) has the //-independent asymptotic behavior e^"^ j ^fa 
at cr ^ CX3 (i.e. z — > 0), the second term in Eq. ([Tl|) is 0(p) and must be kept only when 
= 0. This is the case for the quark distribution at the LO approximation. 
Using the concrete AD values, one has 

/+(z,Q^) = fA, + ^A,Vo(^Lo)e-^+(^)^^° + O(p^o) and 
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9 V ^ 9 
where d+(l) = 1 + 20//(27^o)- 

3.2 the "— " component 

In this case the anomalous dimension is regular and using Eq. (|^)[| one has 

e-An)A,e-^-^-> ^ -el,{\)A,e~'-^'>-o + o(z) 

for PD. 



'^In the Regge-like case (see [|4| and Appendix) fa{x) ^ x ^ ^ z ^ the Eq. (^ has K{\ + X) in its 
r.h.s. with the accuracy 0{z) for an arbitrary A, including the case A = 0. 
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Using the concrete AD values, we have 



4 
9" 

f7{z,Q') = + 0{z), (13) 



f;{z,Q') = --A.e-'^-W^^o + o{z) and 



where rf_(l) = 16//(27/3o). 

Finally we present the full small x asymptotic results for PD and F2 structure function 
at LO of perturbation theory: 



fa{z,Q') = f^{z,Q') + f-{z,Q') and 

F2{z,Q^) = e-f,{z,Q^) (14) 

where f^Jg, fg and /~ are given by Eqs. (|12D and ([T3|) and e = J^i^f/f is the average 
charge square of the / active quarks. 

Let us now describe the main conclusions that follow Eq. (p^). 

• Our LO results for PD coincide with the ones obtained by Mankievitcz et al. in 
Ref. for / = 4 active quarks. The "+" component part coincides with the Ball 



and Forte LO results [0 



The "+" and "— " components are presented explicitly separated. The "— " compo- 
nent ~ Const is negligible at small x (and large Q^) in comparison with Phohi.o'Lo) 
(as it has been already observed in |2^, |12], |T^) and the LO quark distribution 

is "driven" by the influence of the gluons: fq{z,Q'^) ~ {f /9)pfg{z,Q'^) (see also 
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2^, 0, |T5| [13 )• However, at intermediate Q , the "— " component is essential (as 



it was discussed in |TB|, [T^). Thus, in order to give the more general result valid 
for a wide range, we consider PD (see Eq. (|H|)) as the combinations of the "+" 
and "— " components, where every component evolves independently. 

The separation of the singular and regular parts of the AD performed above leads 
to the possibility of avoiding complicated methods for evaluating the inverse Mellin 
convolution or special analyses of DGLAP equations (see |T^ for a review on these 



methods). In our case, we use the exact solution to get the moments of the PD. 
The simple form of the singular part of this exact solution is easily transformed 
to the z-space. The non-singular part is added by the method of replacing Mellin 



convolution by usual product |2J] . In this case the non-singular part in the z-space 



is equal to the corresponding contribution for the first moment n = 1. 

In the following we resume the steps we have followed to reach the small x approximate 
solution of DGLAP shown above: 

• Use the n-space exact solution. 

• Expand the perturbatively calculated parts (AD and coefficient functions) in the 
vicinity of the point n = 1. 
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The singular part with the form 

Aa{n - i)'=e-'^"*^o/("-^) (15) 

leads to Bessel functions in the 2;-space in the form 

^ (fc+i)/2 

^a(-^) h+i[2^/ dsLolnzj (16) 

The regular part B{n) exp {—d{n)sLo) leads to the additional coefficient (see 
and Appendix) 



B{l)exp{-d{l)sLo) + 0{\JdsLo/lnz) 

behind of the Bessel function ([T6| ) in the ^-space. Because the accuracy is 0{\fdsix)jlnz) 
it is necessary to use only the basic term of Eq. (p!6|), i.e. all terms (n — l)'^ in front 
of exp{—d/{n — 1)), with the exception of one with the smaller k value, can be 
neglected. 

• If the singular part at n 1 is absent, i.e. d = in (|T^), the result in the z-space 
is determined by B{l)exp{—d{l)sLo) with accuracy 0{z). 

We would like to stress that the applicability of the above recipe is not limited by 
the order in perturbation theory but by the form of the singular part of the anomalous 
dimensions. At the ffist two orders of perturbation theory the singular part is proportional 
to ~ (n — 1)^^ but this behaviour does not remain at higher orders. The most singular 



terms have been calculated in |2^, ^ . For example, the gluon-gluon AD has the form 



7gg(n,a) = 7(n,a) +0(^«(^— ^ (17) 

where the terms ~ 0(^a{a/{n — l))'^) have been evaluated |^ very recently. 

The BFKL anomalous dimension 7(77,, a) is obtained by solving the implicit equation 

n — 1 ^ ' 

where the characteristic function x(7) has the following expression in terms of the Euler 
\E'-function: 

X(7) = 2M/(l)-M/(7)-vl/(l-7) 
The expansion of 7(n, a) in powers of = 4Ca« gives: 



. + 2.404(^)* + 2.074(^)* + o((;^)') (18) 

which explicitly shows the of the the term ~ (n — 1)"'^ in the fourth order of the expansion. 
Moreover, the dependence ~ (n — 1)~^ has been found |^ when it is considered the third 
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order terms proportional to yO({a/{n — in the r.h.s. of Eq. (p!7|). Note that the 

regular part of this third order coefficient in the perturbative expansion of the AD is only 
partially knownQ although there is the possibility to calculate it completely |34 



Thus, an extension of our recipe beyond the NLO approximation requires the evalua- 
tion of singular terms ~ (n — 1)"'^ {k > 1) which is not easy and out of the scope of this 
work. We restrict ourselves to the first two orders in perturbation theory. 

4 Next-to-leading order 

At the momentum space, the solution of the DGLAP equation has the formQ 

M„(n,Q2) = M+(n,Q2) + M-(n,Q2) and 

M±(n,g2) = M^{n,Q\Ql)exp{-d^{n)s-D±{n)p) _ _ (19) 

= M^{n, Q^, Ql)exp{-{d±s + ]D±p)/{n - l))exp{-d±{n)s - D±{n)p), 

where the new variable p is p = «(Qo) ~" '^(Q^)- 

In comparison with the LO expressions, the NLO AD leads to the following additional 
factors in the moments: 

• The term proportional to p contributes to the evolution part. It can be represented 

by 

D±in) = rf±±(n) - ^rf±(n) (20) 
Po 

and analogously for its singular and regular parts. 



The terms proportional to <y{Q'^) and a(Qo) change the normalization factor 



M^{n,Q\Ql)= (l-dl^{n)a{Q')]M^{n,Ql) + d"^^{n)a{Ql)M^{n,Ql) (21) 



The NLO components of the r.h.s. of Eqs. (|20|, 21) have the form 



d±±{n) = E 4^!'^ 



a,b=q,g 



a=q,g 



®It is known the n=2, 4, 6, 8 and 10 Mellin moments 



^The representation ( [19| ) for the moments M^{n,Q ) is shghtly different from the one in Ref. |^ 
(see Eqs. (2.139)-(2.144) in |Q) because we prefer to separate the "+" and "— " evolutions exponents 
exactly at NLO. The sum Af+(n, Q^) + M^{n, Q^) coincides with the one of pcj]. 
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where rfi^V) = lab{n)/2(3Q and tI^V) are the NLO AD. 

We would hke to stress that the exponent in Eq. (|19]) contains the NLO contribution 
proportional to Eq. (pO|). We have checked this result by direct calculation of NNLO cor- 



rections to the solution of DGLAP equation (following the review of Buras ||30[), arriving 
to terms with the form ~ D'^{n)/2 

The exponential representation of the term proportional to Eq. (pOD is very important 
in the case of the + component because it contains a singular part when n — 1 that could 
make it of the order of the LO contribution, spoiling the perturbative convergence if one 
attempts to expand the NLO part of the exponent. 

Following the steps given in the previous section, one can easily obtain the small x 
behaviour of the PD and F2 at NLO. It has the form: 



fa{z,Q') = f:{z,Q') + f~{z,Q') and 

/-(z,Q2) = A-iQ',Ql)exp{-d^il)s-D4l)p) + 0{z) (23) 

f^{z,Q') = A+{Q',Ql)Io{a)exp{-d^il)s-D+{l)p) + 0{p) (24) 

/+(z,Q2) = A+{Q\Ql)\^{l-dUlMQ'))ph{cr)+20a{Q')Ioia)^ 

■exp{~d4l)s - D4l)p) + 0(p) (25) 

F2{z,Q') = e-(f,{z,Q') + ha{Q')f,{z,Q')) (26) 



where 



a 



AtiQ',Ql) 



2y {dj^s + D^p)lnz , p 

4r 



{d+s + D^p) 
Inz 



1 _ -/a(g2 

/ 
9 



I + 3(1 + -f)aiQt) 



= li^A, + U,^ , A- = A,-20aiQl)A'^ 



a 



2ln{l/z) 



80 



fa{Q' 



(27) 



A 



1 ' 



(2J 



The components of the singular and regular parts of D± have the form: 



++ 



412 



/ 



-20 



d' 







1 (^^ S ^fi^ ^-^-^r 1643 2 r68 13 



rf__(l) 
rf%(l) 



16 / , , 13 r,, 23 13 , 







dL 



;i) = -3(i + ^) . 



(29) 



These corrections are also needed to extend the NNLO analysis of structure functions ||3^ from 
non-singlet to singlet behaviour. 
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The numerical values of these coefficients are listed in Tab. 1 for a different number 
of active quarks. 

The Eqs. (|23|)- (|29|) together with the recipes in the end of the section 3 are the main 
results of this article. 

Looking carefully Eqs. (P3D-(B^), we arrive to the following conclusions: 



Our NLO results coincide with the corresponding of Ball and Forte in Ref. [0 
if one neglects the "— " component, expands our NLO singular terms {p)^Ik+i{o') 
in the vicinity of the point a = aio and ignores the NLO regular terms (i.e. put 
exp{—D+{l)p) = exp{—D_{l)p) = 1 and cancel the terms proportionals to a{Ql) 
into the normalization factors and A'^). We think, however, that this expansion 
is not so correct because it generates NLO corrections of the order of the LO terms. 

The negative sign of the NLO correction in a (see Eq. (pTj)) makes excellent the 
agreement of our result with the parametrization of F2 obtained by De Roeck and De 
Wolf |]36[. Their result is very similar to our LO form of fg in Eq.(p!2|) if one replaces 



slo ^ Slo definition (y) of slo- The value S = 0.708 has been obtained in 

the fit to HI and ZEUS data. Due to 5 < 1, it shows less Q^-dependence than it is 
predicted by perturbative QCD at LO. This slower Q^-dependence may be explained 
naturally by the negative NLO corrections to cr obtained here. 



The behaviour of eqs. (p3|)-(|26|) can mimic a power law shape over a limited region 
of X, Q"^. 



/,(x,Q2)~x-^»''(-'«') and F,{x,Q')r.x-''" 



The quark and gluon effective slopes Yj-^ = — In /^(z, Q^) are reduced by the 
NLO terms that leads to the decreasing of the gluon distribution at small x. For 
the quark case it is not the case, because the normalization factor A^ of the "+" 
component produces an additional contribution undampening as ~ {lnz)~^ . 

The gluon effective slope X^^^ is larger than the quark slope Xg^^ , which is in excellent 

agreement with a recent MRS and GRV analyses [|1^, p!T[| . 

Indeed, because d/dlnx = d/dlnz, the effective slopes have the form. 



Uz,Q^) P /i(a)(l-<_(l)a(Q2)) + 20a(Q2)/o(cr)/p ^ ^ 

fg{z,Q') + {2f)/3a{Q^)-f,{z,Q^) 

The effective slopes X'^-^-^ and X'p2 depend on the magnitudes Aa of the initial PD 
and also on the chosen input values of Ql and A. At quite large values of Q^, where 
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the "— " component is not relevant, the dependence on the magnitudes of the initial 
PD disappear, having in this case for the asymptotic values: 



2\ 



Xeff,as. r,2. ^ /2(a)(l - 4,(l)«(g^)) + 20a (g^)/i (g) /p 
' ^ ' ^ /i(a)(l-4_(l)a(Q2)) + 20a(Q2)/o(a)/p 

\ ^"^'^^ + 111(1/^)' 

where simbol ~ marks approximations obtained by expansions of modified Bessel 
functions In{cr)- These aprroximations should be correct only at very large a values 
(i.e. at very large and/or very small x). 

We would like to note that at LO, where p = plq, the slope X^J.^"''^ [z , Q'^) = 
X'^ff'°-^(^z,Q'^) coincides at very large a with one obtained in (see also [Q) in 

the case of flat input. At the NLO approximation the slope X'p2''^'^{z,Q'^) lies be- 
tween quark and gluon ones but closely to quark slope X'g^^''^^{z,Q^) (see also Fig. 
3). 

Both slopes Xl^^ decrease with decreasing z. A z dependence of the slope should not 
appear for a PD with a Regge type asymptotic {x~^) and precise measurement of 
the slope Xl^'^ may lead to the possibility to verify the type of small x asymptotics 
of parton distributions. 



5 Results of the fits 

With the help of the results obtained in the previous section we have analyzed F2 
HERA data at small x from the HI and ZEUS collaborations separately. Initially 
our solution of the DGLAP equations depends on five parameters, i.e. Ql, xq, Ag, Ag 
and Ay^{nf = 4). In order to keep the analysis as simple as possible we have fixed 
Aj^(nj = 4) = 250 MeV which is a reasonable value extracted from the traditional 



(higher x) experiments and that has also been used by others [|T7|- The initial scale of the 
PD was also fixed into the fits to Qo = 1 GeV^, although later it was released to study 
the sensitivity of the fit to the variation of this parameter. The analyzed data region 
was restricted to x < 0.01 to remain within the kinematical range where our results are 
accurate. Finally, the number of active flavors was fixed to /=4. 

Tab. 2 contains the results of the fits to HI data using Eqs. (|T^) at LO and ([26| ) 



at NLO. The errors used in the calculation of are statistical and systematic added in 
quadrature. 
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Fig. 1 shows F2 calculated from the fit with > 1 GeV^ given in table 1 in comparison 
with HI data. Only the lower bins are shown. One can observe that the NLO result 
lies closer to the data than the LO curve. The lack of agreement between data and lines 
observed at the lowest x and bins suggests that the flat behavior should occur at 
lower than 1 GeV^. 

In order to study this point we have done the analysis considering as a free pa- 
rameter. In exchange, we have fixed Xq to 1. In this case, our formulas must be used 
with caution because during the fit Qg could eventually reach a very small value such that 
the argument of the square root in equation (26) becomes negative when the NLO term 
proportional to p becomes larger than the LO term proportional to s . The kinematic 
region where this problem happens depends on Ql (Qg < Q"^ < Ql), where the upper 
limit Ql also depends on Qg. For example for the sequence Ql = 0.5/0.4/0.3/0.2 GeV^ 

= 0.7/1/2.3/180 GeV^. Thus, in the fits one should be out of this region where more 
terms in the perturbative expansion (NNLO) are needed. 

Comparing the results of the fits in table 3 with those in table 2 one can notice a 
significant reduction in the value of Ag, Qq and the x^- Fig- 1 the better agreement 
with the experiment of the NLO curve is apparent at the lowest kinematical bins. 

Tab. 3 also contains the results of the combined fit to 3 different data set of ZEUS 
1^. The LO and NLO results are compared with data in Fig. 2. The main conclusion is 
that the parameters extracted from fits to HI and ZEUS are very similar. 

The observed decreasing of the gluon density at small agrees with the traditional 
assumption (see |^) that the main contribution at small comes from valence quarks. 
Thus, this non-zero contribution of singlet quark distribution at small seems to mimic 
the valence quark distribution. Indeed, in our approach the nonsinglet quark distribution 
is omitted and this approximation maybe is not so correct at very small Q^. As it 
was shown in the analysis of Ref. based on BFKL dynamics, the nonsinglet quark 
distributions has the following small x asymptotic behavior: 



/32 

xfNs{x,Q^) r^x'-''^^ with a^s-Wyaio, (32) 

where aio is a Q^-independent coupling constant. The numerical estimations with the 
running coupling (see ||3^) lead to the behaviour similar to (|3^ with aLoiQ"^) having 



some effective energy scale A proportional A^^. 

At very small Ql (and x ~ 10~^) the value of aj^s may reach the large value a^s ~ I5 
and the NS quark distributions may start to contribute together with the gluon and the 
singlet distributions. Unfortunately, it is very difficult to add the NS quark distribution to 
our analysis based on DGLAP evolution. Indeed, in agreement with DGLAP dynamics, 
the non-singlet quark distribution with oat^ ~ 1 should be Q^-independent which seems 
contrary to its functional form (^) with the running coupling aioiQ^)- 

Finally with the help of Eq. (pO]) we have estimated the F2 effective slope using 
the value of the parameters extracted from NLO fits to data. For HI data we found 
0.05 < Xp{^ < 0.30 - 0.37 and for ZEUS 0.07 - 0.09 < Xp{^ < 0.31 - 0.34. The lower 
(upper) limits corresponds to = 1.5 GeV^ {Q^ = 400 GeV^). The dispersion in some of 
the limits is due to the x dependence. Fig. 3 shows that the three types of asymptotical 
slopes have similar values. The NLO values of X^p^''^'^ lie between the quark and the 
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gluon ones but closer to the quark slope Yj^'"''^ . These results are in excellent agreement 
with those obtained by others (see references [jl], |1^, |lT], |3^ and also the review Q and 
references therein). 

6 Conclusions 

We have presented the rules to construct the small x form of parton distributions 
having soft initial conditions at the first two orders of perturbation theory. The rules are 
based on the exact n-space solution and lead to the possibihty of avoiding complicated 
methods for evaluating the inverse Mellin convolution or special analyses of DGLAP equa- 
tions (see |]l3l for a review on this methods). We have presented here the PD form at LO 



and NLO approximations of perturbation theory, where the corresponding exact solutions 
in n-space are fully known. Our expressions have quite simple form and reproduce many 
properties of parton distributions at small x, that have been known from global fits. 

We found the very good agreement between our approach based on QCD at NLO ap- 
proximation and HERA data, as it has been observed earlier with other approaches (see 
the review [Q). Thus, the nonperturbative contributions as shadoving effects [^], higher 



twist effects ^y] and others seems to be quite small or seems to canceled between them 
and/or with ln{l/x) terms containing by higher orders of perturbative theory. To clear 
up the correct contributions of nonperturbative dynamics and higher orders containing 
strong ln{l/x) terms, it is necessary as more precise data and futher effors in developing 
of theoretical approaches. 

It is very useful to evaluate the derivatives of F2 and the parton distributions with 
respect to the logarithms of 1/x and directly from eqs.(p3D-(p^). We have presented the 
result for the x logarithmic dependence. The calculation of dF2 / dln^Q"^) and dfa/dln{Q'^) 
is also extremely important in view of the recent claim of disagreement between data and 
predictions in perturbative QCD (see, however, new set GRV analysis where 



this disagreement decreases essentially). We are considering to present this work and also 
predictions for SF in a forthcoming article E3 . 
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7 Appendix 



For reader convenience we present here the illustration]^ of the method to replace the 
convolution of two functions by a simple product at small x. We limite ourselves to the 
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Contrary to Ref . [M we use here the variable z = x/xq. 
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case of regular behavior of kernel moments at n ^ 1. More detailed analysis can be found 
in 



241 and in 42 



Let us to consider the set of PD with have the different forms: 

• Regge-like form /r^z) = z^^ f{z), 

• Logarithmic-like form fhiz) = z~^ln(l/ z)f{z), 

• Bessel-like form fj{z) = z-Hk{2^ dln{l/ z))f{z), 

where f{z) and its derivative f\z) = df{z)/dz are smooth at z = and both are equal 
to zero at z = 1: 

/(I) = /'(I) = 
1. Consider the basic integral with even n > 1: 

Js,in,z) = ® Mz) ^ f'^y'^ h(-) , 

Jz y \yj 

a) Regge-like case. Expanding f{z) near /(O) , we have 



i = R,L,I 



Js,R{n,z) = z' 



dy 
1 



/(0) + -/«(0) + ... + l(^) /«(0) + 



y 



k\ \y, 



n + 5 
1 



f{Q) + 0{z) 



(Al) 



/(O) + 



^ PH0) + ... + i^-/«(0) + 



.n + 5 n + 5 — 1 k\ n + 6 — k 

The second term on the r.h.s. of Eg. ( [AID '^^^ be summed: 

1 



Js,R{n, z) = z~ 



n + S 



f{0) + O{z) 



Tj-jn + Smi + u) ~ 
T{l + u-n-6) 



Because our interest here is limited by the nonsingular case {n > 1), we can neglect 
the second term and obtain: 



Js,R{n,z) = z~ 



f{z) + 0{z'-') 



n + 6 



b) Logarithmic-like case. Using the simple relation z ^hi{l/z) = {d/d6)z ^ we imme- 
diately obtain 



J&An,z) 



z-^ ln(l/z) 
1 



n + 5 



1 - 



n + 5 
1 

n + 5 



n + 5) ln(l/z 
fL{z) + 0{l/ln{l/z)) 



n + 5) ln(l/a; 

f^{z)+0{z'-' 



/(O) + 0{z) 
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c) Bessel-like case. Representing Bessel function in the form 
h{2^dln{l/z)) = J2- 



d— z 



n! r(n + A; + 1) V dS 
and repeating the above analysis, we have 



5=0 



J5,i{ri,z) = 



n + S 



Mz) + 0\ 



d 



2. Consider the integral 



Is(z)^K(z)0f(z)^ f'^k{y) n-) 

■1^ y y 

and define the moments of the kernel K{y) in the following form 

Jo 

In analogy with subsection 1 we have for the Regge-like case: 



dy y'-' K{y) 



^(0) + -^«(0) + ... + l(^) ^«(0) + 
y k\ \y 



1 



where 



Nr,{z)^ f'dy y^-'Kizy) 

J 



The case = l/{n-\-S) corresponds to K{ii) = |/" and has been already considered 
in subsection 1. In the more general cases (for example, = ^(1 + 5) + 7) we can 

represent the "moment" as series of the sort J2m=i ^ / {n + S + m) . 

So, for the initial integral at small x we get the simple equation: 

hAz) = K^+s Rz) + 0{z'-') = K^+s fniz) + 0{z'-') 
Repeating the analysis of the subsections (lb) and (Ic), one easily obtains 



Is,l{'^^ ^) = ^i+s hiz) + O 
h,i{n, z) = fi{x) + O 



ln{l/z) 



and 



d 

\W/z) 



Thus, in the nonsingular case the results do not depend on the shape of the PD but 
the accuracy of the method decreases for a Zn(l/2;) -dependent PD. 
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Tables 



/ 


d+ 


D+ 




^+(1) 


rf_(l) 


D41) 


3 


-4/3 


1180/81 


101/81 


-43.37 


16/81 


1.974 


4 


-36/25 


91096/5625 


61/45 


-45.49 


64/225 


3.108 


5 


-36/23 


84964/4761 


307/207 


-47.73 


80/207 


4.675 


6 


-12/7 


8576/441 


103/63 


-50.05 


32/63 


6.864 



Table 1 . The values of the parameters used in the calculation of the parton distributions 
as a function of the number of flavors 



Q' > 




A, 


•'■q 


x'/ii.o.p. 


LO (HI) 










1 


1.06±0.07 


2.46±0.19 


0.11±0.02 


114/104 


3 


0.96±0.08 


2.53±0.19 


0.12±0.02 


88/92 


5 


0.83±0.09 


2.47±0.18 


0.15±0.02 


40/83 


8.5 


O.BOiO.lO 


2.30±0.18 


0.18±0.03 


23/67 


NLO (HI) 










1 


0.97±0.08 


1.30±0.11 


0.20±0.03 


63/104 


3 


0.91±0.10 


1.31±0.11 


0.22±0.03 


50/92 


5 


0.81±0.10 


1.28±0.11 


0.26±0.04 


27/83 


8.5 


0.84±0.11 


1.22±0.11 


0.28±0.05 


21/67 



Table 2. The result of the LO and NLO fits to HI (1994) data for different low cuts. 
In the fits Ql is fixed to 1 GeV^ 



Aprox. 


A, 




Ql 


X^/n.o.p. 


LO (HI) 
NLO (HI) 


1.10±0.08 
0.83±0.09 


0.35±0.06 
0.21±0.07 


0.55±0.02 
0.55±0.03 


60/104 
45/104 


LO (ZEUS) 

NLO (ZEi;S) 


1.13±0.07 

0.85±0.08 


0.28±0.05 

0.18±0.05 


0.55±0.02 

0.56±0.02 


174/126 
113/126 



Table 3. The results of the fits to HI and ZEUS (1994) data at LO and NLO with 
free. 
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Figure captions 



Figure 1. The structure function F2 as a function of x for different bins. The 
experimental points are from HI The inner error bars are statistic while the outer 
bars represents statistic and systematic errors added in quadrature. The dashed and dot- 
dashed curves are obtained from fits at LO and NLO respectively with fixed Qo = 1 GeV^ 
(see table 2). The solid line is from the fit at NLO giving = 0.55 GeV^ (see table 3). 

Figure 2. The structure function F2. Experimental points are from ZEUS (squares 
from Ref. and diamonds and crosses from two different types of measurements reported 
in Ref. [^). The error bars are displayed as in Fig. 1. Solid (dashed) lines are calculated 
with the parameters given in table 1 from fits at NLO (LO). 

Figure 3. The asymptotical values of effective slopes X^^'^'"''^, X^J^'""^ and X^pl'""^ cal- 
culated at NLO with the parameters from a NLO fit to HI with xq = 1 (see Tab. (3)). 
Lower curves correspond to x = 3 x 10~^ while the upper ones are for x = 10"^. The 
experimental points are from HI IQ. The error bars are displayed as in Fig. 1. 



20 



F2(x,Q') 



F2(x,Q') 



F2(x,Q') 




F2(x,Q') 



F2(x,Q') 



F2(x,Q') 




X 



eff.as 



(Q^) 




